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ABSTRACT 

The  vector  motion  of  severe  tropical  cyclones  (including  storm,  hurri- 
cane/typhoon  stages)  is  forecasted  by  a  numerical  scheme  which  involves 
two  steps: 

a ,  Numerical  geostrophic  steering  of  the  center  of  the  cyclone  using 
the  U.S.  Navy  Fleet  Numerical  Weather  Facility's  (FNWF)  operationally- 
produced  smoothed  isobaric  height  fields,  called  SR.  The  tropical  pertur¬ 
bations  are  steered  in  one-hour  time  steps  up  to  72  hours,  using  winds 
derived  from  the  SR  analysis  dated  closest  to  warning  time.  SR  500  mb.  in 
the  Pacific  and  SR  700  mb.  in  the  Atlantic  gave  the  most  accurate  forecasts 
on  tests  of  ten  northwest  Pacific  typhoons  and  all  five  north  Atlantic  tropical 
storms  and  hurricanes  in  the  period  15  August-1  November  1965.  Forecasts 
were  made  twice  daily,  0600  and  1800  GMT,  during  this  period  using  bert- 
track  information. 

b.  Next,  the  numerical-steering  prediction  is  objectively  modified  to 
adjust  for  bias  (i.e,,  deficiency  in  both  zonal  and  meridonal  motion)  by 
utilizing  errors  made  in  the  most  recent  12-  and  24-hour  numerical-steering 
forecasts .  Several  modes  of  adjustment  are  employed;  the  most  recent 

12- (12-  and  24-)  hour  numerical- steering  bias  yields  the  most  accurate 
correction  of  subsequent  Atlantic  (Pacific)  forecasts  out  to  periods  of  72 
hours.  The  optimal  Naval  Postgraduate  School  (NPGS)  technique  produces 
forecast  errors  ranging  from  an  average  of  4.2  knots  for  12-hour  forecasts 
to  6.2  knots  for  72-hour  forecasts.  The  U.S.  Navy's  official  forecast 
accuracy  is  excelled  by  the  NPGS  scheme  for  all  time  periods. 

Stratification  of  arror  statistics  by  area,  trajectory  and  stage  of  storm. 
Intercomparison  with  ESSA's  NHC-64  technique,  discussion  of  merits  and 


deficiencies  of  the  research  program  relative  to  operational  forecasts  and 
current  experiments  at  FNWF,  are  discussed. 
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A  NUMERICAL  SCHEME  FOR  THE  PREDICTION  OF 
HURRICANE  AND  TYPHOON  MOVEMENT1  • 2 

by 

Robert  J.  Renard 

Department  of  Meteorology  and  Oceanography 
Naval  Postgraduate  School 
Monterey,  California 

1.  INTRODUCTION 

The  Statement  on  Hurricanes  issued  by  the  American  Meteorological  Society 
[l]  indicates  that  the  desired  degree  of  accuracy  in  forecasting  the  position 
of  severe  tropical  cyclones  is  50  miles  or  less  in  a  24-  to  36-hour  period. 

Such  verification  figures  are  far  from  being  realized  at  the  present  time  as  noted 
in  recently  published  error  statistics  from  the  United  States  Fleet  Weather 
Facility,  Jacksonville,  Florida  [10],  the  United  States  Fleet  Weather  Central/ 
Joint  Typhoon  Warning  Center,  Guam  [8],  the  National  Meteorological  Center, 
Environmental  Science  Services  Administration  (ESSA)  [11],  and  ESSA's  National 
Hurricane  Research  Laboratory,  Miami,  Florida  [4,5,6]. 

The  official  hurricane/typhoon  forecast,  based  on  a  careful  consideration 
of  all  the  available  and  pertinent  subjective  and  objective  techniques,  is  be¬ 
coming  increasingly  dependent  on  the  competitive  contributions  from  the 
numerical  approach.  Both  the  United  States  Navy's  and  ESSA's  numerical 
techniques  already  exceed  the  accuracy  of  many  of  the  forecast  schemes 
used  faithfully  by  operational  forecasters  for  many  years  [8,10]  and  yet 
have  potential  for  still  greater  improvement.  Some  of  this  potential  has 
been  realized  recently  by  the  development  of  a  forecast  scheme  using  certain 
numerically-analyzed  operational  products  generated  by  the  United  States 
Navy's  Fleet  Numerical  Weather  Facility  (FNWF),  Monterey,  California. 

When  coupled  with  an  objective  adjustment,  dependent  only  on  the  char¬ 
acteristics  of  the  storm's  recent  trajectory,  the  numerical  scheme  appears 
a  substantial  increase  in  the  accuracy  of  predicting  movement  of 
tropical  cyclones,  as  compared  to  official  forecasts,  for  forecast  intervals 
up  to  72  hours.  The  subject  research  reported  on  here  represents  a 

1 .  Modified  version  presented  at  the  IUGG-AMS  Fifth  Technical  Conference 

and  ^P10®1  Meteorology,  November  20-28,  1967, 
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coordinated  effort  of  the  Naval  Postgraduate  School  (NPGS)  and  FNWF,  Mon¬ 
terey,  California. 

2.  THE  NUMERICAL- STEERING  PROGRAM 

In  addition  to  the  analyses  of  heights  of  mandatory  isobaric  levels ,  FNWF 

operationally  produces  analyses  of  certain  additive  components  of  these 

height  fields  on  a  twice-daily  basis,  0000GMT  and  1200GMT  [3].  This 

unique  numerical  program,  as  developed  for  FNWF  by  Holl  [2],  performs  a 

mathematical  smoothing  of  the  isobaric  height  fields  with  the  degree  of 

smoothing  dependent  on  the  amplitude  and  wave  length  inherent  in  the  iso- 

hyptic  field.  The  arithmetic  difference  between  the  height  field  (Z)  and  the 

smoothed  height  field  (ZgR)  is  called  the  disturbance  field  (zgD)«  Thus,  at 

any  point  on  the  isobaric  surface  Z  «  ZgR  +  ZgD.  The  ZgR  pattern  may  be 

viewed  as  a  space-mean  height  field,  portraying  long  wave  features,  while 

the  Zgp  contours  depict  the  short  or  minor  wave  components  of  the  original 

lsohypses.  Accordingly,  the  Z  field,  void  of  the  disturbance  flow  to  a 

oK 

certain  .degree ,  may  be  used  to  generate  a  current  appropriate  to  the  steering 

of  tropical  cyclones,  which  are  regarded  as  the  disturbance  elements. 

The  foregoing  interrelations  between  Z,  Z_  and  Z  may  be  seen  in 

oK  ou 

figures  la,  b,  and  c,  each  of  which  portrays  a  portion  of  these  isobaric 
height  fields  for  0000GMT,  20  August  1965.  Z  and  ZgR  contours  are  at 
60-meter  intervals  while  Ze_  isolines  are  at  a  30-meter  interval.  Plotted 

DU 

central  values  include  the  units  figure .  Isolines  are  labeled  in  tens  of 

meters;  the  thousands  figure  is  omitted  on  both  the  Z  and  Zc  fields. 

The  nature  of  the  decomposition  analyses  in  the  case  of  severe  tropical 

cyclones  may  be  noted  from  the  situation  just  equatorward  of  Japan.  Typhoon 

Lucy  is  located  at  28.5°N,  140. 2°E  at  map  time  according  to  the  best  track 

position  [7],  Figure  la  shows  that  FNWF's  500-mb  operational  position 

of  Lucy  is  very  close  to  the  best-track  location  with  a  central  500-mb 

height  of  5806  meters.  The  SD  field  emphasizes  the  perturbation  character 

of  Lucy  with  a  minimum  value  of  -75  meters  at  the  typhoon  center.  Thus, 

Z  -  Z  »  Z__  or  5806  *•  (-75)  ■  5881  meters,  which  may  be  verified  from 
SD  SR 

figure  lb. 
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Having  selected  a  steering  field,  namely  the  SR,  geostrophic  winds  are 
computed  to  yield  the  steering  or  basic  current  used  to  forecast  the  motion 
of  the  tropical  cyclone  centers.  Figure  2  is  a  schematic  diagram  illustrating 
a  section  of  the  Zg„  field  appropriate  to  the  tropics  with  superimposed  grid 
points  representative  of  the  FNWF  linear  I,  J  mesh.  The  numerical  compu¬ 
tation  of  the  geostrophic  steering  wind  is  accomplished  by  first  locating  the 
tropical  cyclone  center  to  the  nearest  .1°  latitude  and  longitude.  This 
point  is  identified  as  I,  J  in  figure  2.  Next,  geostrophic  winds  (Vg  J  are 
computed  by  I,  J  components  and  converted  to  latitude  and  longitude  at  each 
of  the  four  points,  I,  J+l;  1-1,  J;  I,  J-l;  and  1+1,  J.  An  average  of  the  four 
geostrophic  winds  is  used  to  steer  the  cyclone  center  for  one  hour.  For  each 
subsequent  hour,  up  to  72  hours,  the  process  is  repeated. 

The  finite  difference  form  of  the  geostrophic  wind  equation  necessitated 
obtaining  height  information  from  a  distance  of  two  mesh  lengths  from  the 
cyclone  center  in  the  cardinal  I,  J  directions  (i.e,  at  I,  J+2;  1-2,  J;  I,  J-2; 
1+2,  J).  One  mesh  length,  D,  (as  I,  J  to  1+1,  J)  is  381  kms  at  60°  latitude, 
which  reduces  to  about  305  km  at  30°  lat. ,  275  km  at  20°  lat. ,  and  240  km 
at  10°  lat. 

A  potential  problem  with  the  Coriolis  parameter,  used  in  geostrophic 
wind  computations  at  low  latitudes,  was  avoided  by  using  a  modified  form 
of  the  sin  function  for  latitudes  less  than  30°  lat.: 

Mod  sin  8=  2[(.25  sin  8  +  .25)^  +  .25  sin  8)]  . 


The  function  is  graphed  in  figure  3.  The  magnitude  of  mod  sin  8  ranges  from 
.125  at  0°  lat.  to  .53  at  30°  lat.  The  lower  limit,  .125,  is  the  value  of 
the  sin  8  at  7.2°.  Along  with  using  an  average  geostrophic  wind,  as  de¬ 
scribed  above,  the  adjustment  of  the  Coriolis  parameter  may  be  viewed  as 
a  further  reduction  of  the  steering  wind  relative  to  the  true  value  at  the 
position  of  the  tropical  cyclone  center. 

The  steering  section  of  the  forecast  program  was  written  for  operation 
on  the  Control  Data  Corporation's  1604  digital  computer.  Both  NPGS  and 
FNWF  computers  were  utilized  for  these  computations. 
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3.  COMPOSITION  OF  THE  TEST 

Due  to  limited  manpower  and  computer  time ,  only  those  named  North 
Pacific  and  Ncrth  Atlantic  tropical  cyclones  in  existence  during  the  period 
15  August  to  1  November  1965  were  incorporated  into  the  test.  The  sample 
included  all  1965  Atlantic  hurricanes  and  tropical  storms  (Anna,  Betsy, 

Carol,  Debbie,  Elena)  and  ten  Pacific  typhoons  (Lucy,  Mary,  Olive,  Rose, 
Shirley,  Trix,  Virginia,  Bess,  Carmen  and  Della).  Only  0600  and  1800  GMT 
best  track  positions,  as  given  in  [7]  and  [9]  were  used,  except  for  hurri¬ 
cane  Carol,  in  which  case  0000  and  1200  GMT  positions  were  employed. 
Cyclone- position  forecasts,  made  in  one-hour  time  steps,  were  printed  out 
for  each  12-hour  forecast  interval  up  to  48  hours  and  at  72  hours.  De¬ 
pression,  tropical  storm,  hurricane/typhoon  and  extrairopical  stages  were 
included  if  the  position  was  listed  in  the  annual  summaries .  For  the  Atlantic 
area,  79%  of  the  forecasts  were  made  during  storm  and  hurricane  stages, 

13%  from  depression  and  8%  from  extratropical  stages.  In  the  Pacific  the 
vast  majority,  98%,  were  from  the  storm/typhoon  stages;  the  remainder, 

2%,  were  depressions. 

It  is  important  to  note  that  best-track  cyclone  positions  were  used  in 
generating  the  forecasts  up  to  72  hours  using  geostrophic  steering  winds 
computed  from  a  single  SR  analysis  dated  at  best-track  time  plus  six  hours . 
Thus ,  in  order  to  forecast  the  movement  of  a  cyclone  positioned  at  I ,  J 
(figure  2)  at  0600  GMT  (1800  GMT)  the  analyzed  SR  field  for  1200  GMT 
(0000  GMT)  was  used.  However,  in  the  case  of  Carol,  initial  hurricane 
position  and  SR  steering  flow  were  for  the  same  synoptic  time. 

Although  a  different  combination  of  best-track  time  and  time  of  num¬ 
erical  analysis  may  have  been  more  operationally  realistic  the  major  effort 
to  this  point  was  directed  toward  establishing  the  feasibility  of  using  the 
SR  field  to  derive  a  steering  current.  In  addition,  initial  efforts  were 
concentrated  on  the  12-hour  forecasts  for  which  time-mean  SR  steering 
winds  are  appropriate.  From  an  operational  point  of  view  the  SR  analyses 
used  in  the  test  may  be  regarded  as  "perfect"  12-hour  SR  prognostic  fields 
initiated  six  hours  before  warning  time.  Using  SR  analyses  rather  than 
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prognoses  to  derive  the  steering  current  has  the  advantage  that  the  forecast 
errors  do  not  include  contributions  from  the  deficiencies  of  a  numerical 
prognostic  model . 

4.  PRELIMINARY  FINDINGS 

Some  of  the  preliminary  findings  which  established  the  format  for  uni¬ 
formly  processing  the  forecast  data  to  72  hours  is  important  and  will  be  out¬ 
lined  here.  The  initial  experiments  were  designed  to  provide  answers  to  the 
following  questions: 

a)  Which  SR  field (s)  will  generate  the  most  accurate  forecasts  con¬ 
sidering  both  Atlantic  and  Pacific  areas,  and 

b)  does  the  forecast  accuracy  deteriorate  as  the  forecast  interval  ls 
increased  to  72  hours  ? 

Table  1  indicates  the  average  12-hour  forecast  errors  (in  nautical  miles) 
for  each  cyclone  which  resulted  from  applying  the  numerical  steering  pro¬ 
gram  (section  2)  to  SR  fields  from  several  selected  isobaric  levels  (1000, 
700,  500,  200  mb)  and  layers  (1000/500,  1000/200,  500/200  mb).  The 
official  forecast  errors  from  [7]  and  [9]  are  also  shown.  In  the  case  of  the 
Pacific  cyclones  the  overall-average  official  forecast  error  (70n.mi.)  was 
derived  from  a  linear  extrapolation  of  published  error  data  at  24,  48  and  72 
hours.  The  number  of  forecasts  made  are  shown  in  parentheses.  Although 
a  forecast  for  every  0600  and  1800  GMT  best-track  position  given  in  [7] 
and  [9]  was  attempted,  missing  or  unretrievable  SR  data  disallowed  some 
cases,  especially  those  involving  200*mb  data.  Superscripts  "1"  and  "2” 
indicate  the  relative  merits  of  the  two  SR  fields  yielding  the  most  accurate 
12-hour  forecasts,  for  each  named  storm  and  for  the  overall  average. 

Following  are  pertinent  conclusions  to  be  drawn  from  table  1 . 

a.  Considering  both  oceans  collectively,  SJl  500  performed  best 
(i.e.  least  forecast  error),  SR  700  second  best  and  SR  1000/500  a  close 
third. 

b.  In  the  Atlantic,  on  the  average,  SR  700  is  best,  followed  by 

SR  1000  and  SR  1000/500.  Tha  majority  of  individual  cyclones  (three  out  of 
five)  behaved  like  the  overall  average . 
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c.  In  the  Pacific,  on  the  average,  SR  500  is  best,  followed  by  SR  700 
and  SR  1000/500,  each  with  similar  results.  The  vast  majority  of  storms 
(seven  out  of  ten)  behaved  like  the  overall  average. 

d.  SRs  at  the  higher  levels  and  layers  (SR  200,  SR  1000/200  and  SR 
.500/200)  generally  yielded  the  poorest  results. 

e.  Even  considering  the  best  numerical  steering  result  for  the  Atlantic 
(SR  700)  and  Pacific  (SR  500),  the  official  forecast  accuracy  is  superior  by 
more  than  10%. 

In  view  of  the  good  performance  of  SR  500  in  both  oceans  this  field 
was  selected  to  test  the  feasibility  of  SR  steering  computations  for  forecast 
intervals  out  to  72  hours.  Table  2  shows  the  result.  Format  of  error  data 
for  this  and  most  of  the  following  tables  is  similar  to  table  1 .  Official 
36-hour  errors  (from  [7]  and  [9])  are  linear  extrapolations  from  published 
enror  data  at  24,  48  and  72  hours  in  the  Pacific  and  12,  24,  48  and  72  hours 
In  the  Atlantic.  Table  3  is  offered  as  an  aid  in  summarizing  the  results 
shown  in  table  2 . 

‘  s 

a.  For  each  of  the  forecast-time  intervals  tested,  official  results 
excel  the  NPGS  steering  system,  except  for  72  hours  in  the  Atlantic. 

b.  The  forecast  error,  in  nautical  miles/hour  of  forecast  interval, 
generally  decreased  or  held  steady  with  time  for  the  NPGS  system,  while 
the  official  error  figure  increased  or  held  steady  with  time. 

c.  In  a  relative  sense,  the  NPGS  scheme  shows  improvement  compared 
to  the  official  out  to  72  hours,  especially  in  the  Atlantic.  This  may  be 
seen  from  the  calculations  of  NPGS-OFFICIAL  errors  at  each  forecast  in¬ 
terval  . 

A  consideration  of  the  results  shown  in  tables  1-3  suggested  continued 
and  more  extensive  experimentation  with  the  steering  technique.  However, 
before  embarking  on  further  testing  and  possible  modification  of  the 
numerical  scheme  for  periods  out  to  72  hours,  some  further,  but  limited, 
checks  on  the  apparent  merits  of  SR  700  in  the  Atlantic  and  SR  500  in  the 
Pacific  were  attempted. 
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Due  to  the  favorable  performance  of  both  SR  500  and  SR  700,  consider¬ 
ing  both  oceans,  and  similarly  good  behavior  of  these  levels  from  other  but 
related  forecast  techniques  [6],  the  two  levels  were  intercom  pared  at  other 
than  the  12-hour  forecast  intervals . 

Table  4  shows  the  relative  merits  of  SR  700  and  SR  500  for  Atlantic  fore¬ 
casts  up  to  a  period  of  72  hours.  Error  statistics  are  shown  for  all  forecast 
data  as  well  as  for  a  homogeneous  sample  of  forecast  times.  For  example, 

93  SR  700  and  89  SR  500  forecasts  were  possible  for  the  12-hour  interval 
while  a  maximum  of  87  forecast  times  were  common  to  the  SR  500  and  SR  700 
computations.  Results  indicate  the  excellence  of  the  SR  700  forecasts  out 
to  48  hours  and  SR  500  thereafter.  The  official  forecast  is  superior  to  the 
optimal  NPGS  system  in  the  Atlantic  (i.e.  SR  700)  at  12  and  24  hours,  while 
the  official  is  worse  than  the  optimal  NPGS  scheme  at  72  hours.  At  36  and 
48  hours  official  and  NPGS  accuracy  are  nearly  equivalent. 

Testing  of  the  relative  worth  of  SR  700  and  SR  500  in  the  Pacific  was 
limited  to  12  and  24  hours  only  (table  5).  SR  500  maintained  a  lead  in 
accuracy  over  SR  700  through  24  hours  but,  as  in  the  Atlantic,  official  fore¬ 
cast  accuracy  surpassed  the  NPGS  SR  500  geostrophic  steering  forecasts. 

The  preliminary  findings  displayed  in  tables  1-5  led  to  the  decision  to 
use  SR  500  in  the  Pacific  and  SR  700  in  the  Atlantic,  exclusively,  for  all 
subsequent  testing  of  the  numerical  tropical  cyclone  steering  scheme. 

5 .  THE  MODIFIED  NPGS  FORECAST  SCHEME  USING  A  CORRECTION 
FOR  BIAS  IN  THE  NUMERICAL  STEERING  COMPUTATION 

A  typical  example  of  using  the  SR  700  geostrophic  steering  computation 
is  shown  in  figure  4,  The  best  track  positions  (from  [9])  and  the  24-hour 
forecast  positions,  using  numerical  steering  only,  denoted  by  F ^  ,  are 
plotted  for  hurricane  Elena.  Numbered  positions  on  the  best  track  are  at 
12-hour  intervals,  starting  with  1800  GMT  12  October  1965,  which  is 
labeled  position  "O".  Corresponding  numbers  along  the  F forecast  track 
refer  to  the  same  times  as  those  on  the  best  track.  For  example,  "5"  on  the 
best  track  is  the  cyclone-center  position  for  0600  GMT  15  October  1965 
while  position  "5"  along  F ^  is  the  forecast  cyclone-center  position  for  the 
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same  time.  24-hour  forecast  position  "5"  was  generated  from  best  track 
position  "3"  at  0600  GMT  14  October  1965,  using  the  SR  700  analysis  at 
1200  GMT  14  October  1965. 


Figure  4  indicates  that  the  forecast  and  best  tracks  are  similar  in 

i 

shape  but  positions  at  given  times  are  not  identical.  This  feature,  common 
to  most  hurricanes  and  typhoons  considered  in  this  research,  may  be  de¬ 
scribed  as  a  consistent  deficiency  in  both  zonal  and  meridional  components 
of  the  numerical-steering  forecast.  Thus,  the  vector  error  between  forecast 
and  best-track  positions  represents  a  bias  which  may  be  used,  with  ad¬ 
vantage,  as  a  correction  or  modification  to  the  subsequent  numerical- steeuog 
forecasts.  \ 

Figure  5  schematically  indicates  the  mechanics  of  applying  one  type  of 
correction  for  the  bias  in  the  numerical  steering  forecast  (F  ).  For  this  fig- 

XX 

ure  and  in  the  discussion  which  follows  the  subscripts  associated  with  T  and 

F  (as,  -12,  0,  12,  24,  etc.)  refer  to  time  before  or  after  "0",  where  "0"  is 

the  time  at  which  the  forecast  in  question  is  made  (i.e,  warning  time).  The 

solid  lines  connect  best-track  positions  for  -24,  -12,  0  and  24  hours,  while 

F  and  F_,  indicate  the  24-hour  numerical  forecast  positions.  E„  .  is  the 
O  24  24 

vector  error  of  the  24-hour  forecast  made  from  the  T  ^  position.  This 

forecast  error,  known  at  time  "0",  is  then  employed  as  a  correction  to  the 

24-hour  forecast  made  at  time  "0" .  Such  a  procedure  generates  modified 

24 

numerical  forecast  positions,  as  in  figure  5.  The  superscripts  on  F, 
as  used  in  figure  5,  and  in  the  text,  figures  and  tables  which  follow  refer 
to  the  forecast  interval  from  which  the  correction  for  numerical-forecast  bias 
was  selected.  Hence,  superscript  "24"  refers  to  use  of  the  most  recent 
24-hour  numerical  forecast  error  as  a  correction  to  the  numerical  steering 
forecast  made  at  time  "0"  .  Applied  to  a  48-hour  forecast,  the  scheme 

48 

symbolized  in  figure  5  yields  a  modified  forecast  position  designated  F^g  . 

It  is  to  be  noted  that  as  the  forecast  interval  increases  so  does  the 
necessary  time  lag  for  application  increase.  Thus,  in  order  to  make  a 
modified  72-hour  forecast  one  72-hour  forecast  period  must  pass  before 
If  is  known.  This  limits  quite  severely  the  totality  of  application  in  the 

r  m 

forecasting  of  tropical  cyclones. 


V  V 
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In  view  of  the  difficulty  just  mentioned  a  scheme  similar  to  that  shown 

in  figure  5  was  developed,  but  in  this  case  only  the  most  recent  12-hour 

numerical  forecast  errors  were  used,  regardless  of  forecast  interval.  Figure 

6  shows  a  modified  24-hour  forecast  made  at  time  "0"  employing  numerical 

steering,  F_.,  plus  a  correction  for  bias,  2E  _.  The  modified-forecast 

position  is  designated  F,  .  .  When  applied  to  a  48-hour  (72-hour)  forecast 

“  •  12  12  ^  ^ 

the  modified-forecast  position  is  F^g  (F72)  and  4E12(6El2)  is  the  appropriate 
vector  correction  for  bias. 

The  two  modes  of  modifications  just  outlined  are  applied  to  24-hour 

forecasts  of  Elena  and  are  shown  in  figure  7.  It  is  obvious  that  the  modified 

12 

forecasts  result  in  cyclone-position  forecasts  superior  to  F24  with  the  F^ 
scheme  best. 

For  forecast  intervals  beyond  12  and  24  hours  two  other  modification 
schemes,  involving  the  bias  correction,  were  used.  Figures  8  and  9  in¬ 
dicate  these  modes  as  applied  to  36-hour  forecasts  made  at  time  "0”.  The 
scheme  in  figure  8  is  exactly  lik'*  that  in  figure  6  only  24-  rather  than  12- 
hour  values  of  E  are  employed.  Figure  9  shows  a  scheme  for  which  the 
most  recent  12-  and  24-hour  numerical  forecast  errors  are  used,  the  former 
weighted  twice  that  of  the  latter. 

Examples  of  the  type  of  modifications  portrayed  in  figures  8  and  9  are  . 

shown  in  figure  10  for  48-hour  forecasts  of  typhoon  Carmen.  Though  not 

24 

entirely  obvious  at  this  point  the  F^g  system  gives  the  optimal  forecast 

i 

track  for  Carmen. 

Each  one  of  the  adjustment  schemes  described  in  figures  5,  6,  8  and  9 
was  applied  to  each  Pacific  typhoon  for  every  possible  0600  and  1800  GMT 
forecast  time  while  application  to  Atlantic  tropical  storms  and  hurricanes 
was  somewhat  more  limited.  Tables  6  to  16  give  the  basic  results  using 
SR  500  in  the  Pacific  and  SR  700  in  the  Atlantic.  Wherever  possible, 
averages  for  homogeneous  sets  of  forecast  data  are  shown.  In  the  case  of 
official  forecasts,  such  comparisons  are  limited  due  to  different  forecast 
times  in  the  Atlantic  (04,  10,  16,  22  GMT)  and  the  availability  of 
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individual  official  forecasts  for  24  and  48  hours  only  in  the  pacific.  Various 
aspects  of  these  error  statistics  are  summarized  below: 

12 

12  hours:  Atlantic  (table  6):  The  modified  forecasts  (F12)  represent  a 

48%  improvement  over  the  forecasts  as  well  as  significantly  excelling 

the  accuracy  of  the  official  forecasts.  This  is  true  for  the  overall  average 

and  for  each  individual  storm. 

Pacific  (table  7):  The  situation  in  the  Pacific  is  similar  to 
12 

the  Atlantic  with  F^  representing  a  34%  improvement  relative  to  F^2  and  a 
23%  increase  in  accuracy  over  the  estimated  official-forecast  score. 

24  hows:  Atlantic  (table  8):  The  12-hour  numerical-steering  forecast 

/ 

bias  is'most  significant  for  the  24-hour  forecast,  a  not  unexpected  result 

12 

irKview  of  the  relation  of  F12  to  F12 .  Again  the  official-forecast  error  is 

considerably  greater  than  that  of  the  NPGS  optimal  scheme,  namely  F*jj. 

Pacific  (table  9):  Unlike  the  Atlantic,  application  of  the 

24 

24-hour  bias  correction  yields  the  optimal  scheme,  F24,  with  the  official 

forecast  error,  150  n.mi. ,  again  in  considerable  excess  of  the  108  n.mi. 

error  for  a  homogeneous  sample  of  61  F^  forecasts.  Every  cyclone,  except 

Rose,  shows  the  F24  error  less  than  the  official  figure. 

36  hours:  Atlantic  (table  10):  Although  ail  three  schemes  equal  or 

12 

surpass  the  extrapolated  official  forecast  accuracy,  Foc  is  optimum,  again 

O  0 

emphasizing  the  importance  of  the  most  recent  history.  Anna  provides  a 

minor  exception  to  this  trend . 

Pacific  (table  11):  The  importance  of  the  24-hour  forecast 

24 

to  the  modified  36-hour  forecasts  is  seen  from  the  performance  of  Foc  with 

o  b 

an  average  error  of  170  n.mi.  (63  cases)  compared  to  225  n.mi.  for  the 
official  forecasts. 

48  hours:  Atlantic  (table  12):  At  this  point  in  the  range  of  forecast 

intervals  considered,  a  correction  x.r  bias  taken  from  the  same  interval  as 

48 

the  forecast  (i.e.  48  hours)  is  detrimental.  That  is,  the  errors  using  F._ 

12  40 
are  greater  than  for  F.„.  The  F  scheme  continues  to  be  best  while  the 

48  xx 

official  and  optimal  NPGS  schemes  are  producing  results  quantitatively 
more  similar  than  for  shorter  forecast  periods. 
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Pacific  (table  13):  The  Pacific  sample  continues  to  behave 

12  24 

differently  than  the  Atlantic  with  the  most  complex  modification,  F^g'  , 
yielding  the  best  forecast  accuracy,  although  all  four  types  of  bias 
corrections  improve  upon  F^g.  The  official  forecast  continues  to  be  ex¬ 
celled  by  the  optimal  NPGS  scheme, 

72  hours:  Atlantic  (table  14):  Comments  for  the  Atlantic  at  48-hours 
are  true  for  the  72-hour  forecast  interval  as  well  except  that  bias  corrections 
of  any  type  considered  did  not  improve  upon  F72»  However,-  official-forecast 
accuracy  is  still  surpassed,  namely  by  F^* 

Pacific  (table  15):  This  is  the  most  difficult  table  to  interpret 

since  the  averages  for  non-homogeneous  sets  of  forecast  data  suggest 

12  24  24 

F^  Is  best  while  the  homogeneous  sets  of  data  indicates  F^  is  optimum. 

But,  since  the  number  of  cases  is  relatively  small  for  the  homogeneous  sets, 
the  latter  figures  cannot  be  regarded  as  significant. 

Summarizing  information  for  tables  6-15  is  shown  in  table  16.  The 
forecast  error  per  unit  of  time,  using  the  cjrtimal  NPGS  scheme,  increases 
with  time,  particularly  so  in  the  Atlantic.  However,  the  NPGS  system  always 
surpasses  the  official  forecast  accuracy  although  the  ratio  generally  decreases 
with  increasing  time.  In  addition,  table  16  gives  information  on  the  dis¬ 
tribution  of  NPGS  forecast  errors,  using  the  optimal  scheme.  In  the  Atlantic, 
on  the  average,  about  2/3  (1/2)  of  the  forecast  errors  lie  within  3  kts.  of 
the  average  forecast  error  through  36  (for  48  and  72)  hours .  The  dispersion 
of  errors  is  considerably  less  in  the  Pacific  where  3/4  (2/3)  represents  the 
corresponding  number  of  cases  for  12,  24  and  36  (48  and  72)  hours.  Con¬ 
sidering  both  oceans,  the  remaining  1/3  of  the  cases  are  about  evenly 
distributed  between  the  very  large  (greater  than  average  plus  3  knots)  and 
very  small  (less  than  average  minus  3  knots)  forecast  errors.  It  is  also 
evident  from  the  listing  of  optimal  schemes  that  the  short-term  peculiarities 
(i.e.  12,  24  hours)  in  cyclone  trajectories  have  long-term  application  (up 
to  72  hours)  in  the  modified  forecast  procedure. 

Figures  11  and  12  show  Atlantic  and  Pacific  examples  of  forecasts  to 
72  hours,  each  made  from  a  given  synoptic  time.  The  best  track,  the  NPGS 

11 


numerical-steering  forecast  track  (F  )  and  the  optimal  NPGS  forecast  track 

.  XX 

(F  p  )  are  shown  for  the  two  cases  portrayed.  Additionally,  the  available 
official-forecast  positions  are  indicated.  The  inadequacy  of  the  numerical- 
steering  forecast  relative  to  the  modified  forecast  is  clearly  indicated  in 
both  the  Debbie  and  Rose  figures.  The  reasonable  continuity  of  successive 
forecast  positions,  12  to  72  hours,  using  the  optimal  scheme  is  evident. 

The  extreme  disparity  which  may  occur  between  the  official  and  the  modified 
numerical  schemes  is  also  shown  in  the  case  of  Debbie. 

6.  NHC-64  vs  NPGS  OPTIMAL  FORECAST  SCHEME 

A  further  evaluation  of  the  NPGS  forecast  errors  was  made  through  an 
intercomparison  with  the  NHC-64  statistical  technique  [4,5,6]  as  developed 
by  the  National  Hurricane  Research  Laboratory,  Miami,  Florida.  Table  17 
shows  results  for  12-,  24-,  36-  and  48-hour  forecasts.  Since  the  NHC-64 
forecasts  were  made  at  0000  and  1200  GMT,  an  average  of  the  errors  from 
0600  and  1800  GMT  optimal  NPGS  forecasts  were  compared  to  each  NHC-64 
forecast  considered.  This  is  the  closest  approach  to  homogeneity  that 
could  be  made  here.  Carol  is  an  exception,  since  0000  and  1200  GMT 
NPGS  forecasts  were  computed  making  this  storm's  sample  truly  homogeneous 
with  the  NHC-64  cases.  Average  errors  indicate  the  NPGS  optimal  scheme 
excelled  NHC-64  for  12  and  24  hours  while  the  latter  surpassed  the  former 
at  48  hours,  although  even  here  two  out  of  three  storms  representing  50% 
of  the  forecasts  favor  the  NPGS  scheme.  36-hour  statistics  yield  incon- 
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elusive  results.  72-hour  forecasts  were  not  available  from  Miami.  For 
Carol,  whose  forecast  times  exactly  matched  those  of  NHC-64,  the  results 
show  the  NPGS  scheme  to  be  an  improvement  over  NHC-64  at  all  forecast 
intervals . 

7.  STRATIFICATION  OF  ERROR  STATISTICS  BY  AREA,  TRACK  AND  STORM  STAGE 

Atlantic:  The  Atlantic  area  was  divided  into  three  zones.  A,  B,  and  C, 
in  accordance  with  a  similar  division  used  by  the  NHC  group  at  Miami  [6], 
See  figure  13.  Area  A  represents  the  Atlantic  area  generally  east  of  60°W,; 
while  B  covers  the  Caribbean,  Gulf  of  Mexico  and  Atlantic  areas  south  of 
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30°N,and  north  and  west  of  60°W.  C  encompasses  the  eastern  United  States 
and  ocean  areas  immediately  to  the  east  which  are  north  of  30°N. 

Table  18  shows  the  predominance  of  cyclone  positions  in  areas  A  and  B 
and  the  superiority  of  optimal  NPGS  forecast  accuracy  in  the  latter  compared 
to  A  and  C  for  all  forecast  time  intervals .  In  the  case  of  72-hour  forecasts, 
errors  in  areas  A  and  C,  collectively,  average  about  100%  greater  than  those 
in  area  B.  Such  statistics  compare  well  with  findings  by  Tracy  [6]  using 
the  NHC-64  technique. 

The  distribution  of  forecast  errors  relative  to  path  is  also  quite  inter¬ 
esting.  Without  exception,  forecast  errors  are  less  for  all  forecast  inter¬ 
vals  for  cyclone  stages  before  the  time  of  recurvature.  After-recurvature 

in 

areas  for  the  storm/hurricane  stages  are  most  frequentlyAC  and  northern 
sections  of  A  as  shown  in  [9]. 

The  interrelationship  of  area  and  path  are  also  manifest  in  the  error 
statistics  relative  to  cyclone  stage.  Table  18  indicates  that,  collectively, 
the  intensifying  tropical  depression  (TD  dev)  and  tropical  storm  (TS)  are 
associated  with  the  most  reliable  results  when  using  the  optimal  NPGS 
scheme.  These  stages  are  generally  before  recurvature  and  in  area  B  or 
southwestern  A.  Hurricane  (H)  statistics  are  next  best,  partly  due  to  in¬ 
clusion  of  some  after-recurvature  cases  in  areas  A  and  C.  The  extra- 
tropical  (EXT)  and  dissipating  tropical  depression  TD(dis)  stages  should  be 
combined  as  case  histories  in  [9]  lead  to  the  conclusion  that  the  differ¬ 
ences  between  the  two  stages  are  quite  tenuous.  All  in  ail  EXT  and  TD(dis) 
categories  perform  poorest  and  represent  after-recurvature  cases  in  area  A 
for  the  most  part. 

Pacific:  Table  19  for  Pacific  typhoons  shows  the  breakdown  by  path 
only.  Area  analysis  has  not  received  the  same  focus  as  in  the  Atlantic  and 
analysis  by  stage  from  published  1965  storm  data  [7]  does  not  discriminate 
sufficiently  between  cyclone  categories  to  warrant  analysis  like  that  pre¬ 
sented  for  the  Atlantic.  Before  recurvature,  error  values  are  much  less  than 
the  overall  average  official  errors  while  after  recurvature  the  optimal  NPGS 


errors  jump  by  as  much  as  100%.  These  results  closely  parallel  the  Atlantic. 

8.  CONCLUSIONS 

The  NPGS  scheme  for  forecasting  tracks  of  tropical  storms,  hurricanes 
and  typhoons  is  objective,  numerical,  easy  to  apply  and  readily  adaptable 
to  field  use.  The  errors  for  forecast  intervals  up  to  72  hours  are  consistently 
below  those  from  most  other  well  known  subjective  and  objective  techniques. 

Part  of  the  success  of  the  NPGS  scheme  in  relation  to  the  official  and 
NHC-64  forecasts  may  be  ascribed  to  the  following. 

a.  Best-track  vice  operaticnal-track  positions  were  used  as  initial 
cyclone  locations  from  which  NPGS  forecasts  were  generated,  while  the 
operational  positions  are  germane  to  the  official  and  NHC-64  statistics. 
However,  all  three  techniques  used  best-track  data  for  verifications. 

Thirteen  miles  is  the  average  difference  between  aircraft  reconnaissance 
and  best- track  locations  in  the  Pacific  in  1965  [7],  Such  a  difference  rep¬ 
resents  a  range  from  about  25%  to  4%  of  the  magnitudes  of  the  forecast 
errors  for  periods  from  12  to  72  hours,  respectively.  This  factor  does  not 
appear  to  change  the  conclusions  cited  to  this  point. 

b.  Perhaps  more  serious  than  (a)  above  is  the  following.  In  the  case 
of  JTWC/FWC  Guam  the  operational  positions  at  forecast  times  (0600  and 
1800  GMT)  are  determined  by  three-  to  twelve-hour  forecasts  from  fixes 
determined  by  recent  land  radar  and/or  aircraft  reconnaissance  observations 
or  by  surface  upper  air  analyses.  Such  a  procedure  puts  the  official  fore¬ 
cast  at  a  disadvantage  compared  to  the  research  program  used  here.  The 
magnitude  of  the  disadvantage  is  difficult  to  assess. 

c.  As  noted  in  Section  3,  SR  analyses,  six  hours  after  initial  time, 
were  used  to  forecast  cyclone  tracks  out  to  72  hours.  This  is  not  operation¬ 
ally  realistic  and  may  have  contributed  somewhat  to  the  success  of  the  NPGS 
scheme,  particularly  in  the  short-period  forecasts  as  12  and  24  hours. 

Balancing  the  scale  in  favor  of  the  relative  merits  of  the  NPGS  scheme 
is  the  recent  operational  experience  of  JTWC/FWC  Guam.  In  the  summer 
of  1967  FNWF  began  an  experimental  numerical  tropical  cyclone  steering 
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program  which  utilizes  the  SR  fields  in  essentially  the  same  way  as  the  re¬ 
search  program  outlined  here.  The  movement  forecasts  are  produced  separate¬ 
ly  from  SR  analyses  and  prognostic  fields.  Guam  has  used  these  numerical 
steering  forecasts  along  with  corrections  for  bias  in  the  manner  just  described. 
Preliminary  indications  suggest  that  the  accuracy  of  24-hour  forecasts, 
accomplished  under  operational  real-time  conditions  is  commensurate  with 
that  shown  in  this  paper,  as  performed  under  a  post-season  research 
environment.^  Definite  statements  on  this  matter  await  extensive  post¬ 
season  analysis . 

9 .  FINAL  REMARKS  AND  AVENUES  FOR  FURTHER  RESEARCH 
The  merits  of  the  bias  correction  are  derived  from  the  information  content 
inherent  in  the  recent  behavior  of  the  storm  relative  to  the  numerical  scheme 
used  to  predict  it.  This  is  a  simple,  however  unique,  application  of  con¬ 
tinuity.  As  such,  the  correction  for  bias  using  SR  analyses  only  may  be 
viewed  as  serving  one  or  more  of  the  following  purposes .  It  compensates 
for  (a)  the  use  of  an  improper  steering  field  and/or  derived  current,  and/or  (b) 
the  use  of  an  inappropriate  level  or  layer  in  the  SR  steering  field  and/or 
(c)  erroneous  information  in  the  particular  SR  field  selected  as  the  steering 
medium,  and/or  (d)  changes  with  time  in  the  SR  steering  field.  The  last 
point  is  tantamount  to  stating  that  the  correction  for  bias,  especially  at 
increasing  forecast  intervals,  substitutes  for  movement  and  development  in 
e  SR  steering  current,  but,  of  course,  with  lag.  Since  prognostic  fields 
are  imperfect,  especially  in  the  tropics,  the  procedure  of  using  a  bias 
correction  to  approximate  changes  in  the  SR  field  may  be  preferable .  Ex¬ 
periments  are  being  conducted  at  both  the  FNWF  and  the  NPGS  to  determine 
the  merits  and  deficiencies  of  using  SR  analyses  only  or  SR  analyses  and 
prognoses  in  combination,  to  generate  forecasts  of  tropical-cyclone  move¬ 
ment.  Perhaps  the  temporal  deterioration  of  the  information  content  in  the 
initial  SR  analyses  suggests  using  a  relatively  reliable  short  period  SR 
prognostic  field,  as  the  36-hour,  for  cyclone  forecasts  from  36  to  72  hours. 

3  Private  communication  with  personnel  at  JTWC/FWC  Guam 
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More  directly,  a  consistent  bias  in  the  numerical  steering  program  strong¬ 
ly  suggests  tuning  the  steering  field  or  its  derived  current  to  the  movement 
of  tropical  cyclones.  In  other  words ,  changes  may  be  made  to  the  mathe¬ 
matical  smoothing  program  to  allow  increased  meridional  steering  components 
as  well  as  magnification  of  the  basic  zonal  current. 

Further,  utilizing  the  geostrophic  SR  wind  at  the  point  of  the  storm  cen¬ 
ter  instead  of  a  mean  geostrophic  wind  from  the  area  surrounding  the  storm  is 
likely  to  give  some  increase  in  the  steering  values.  Such  a  modification  is 
already  a  part  of  the  present  FNWF  experimental  tropical-cyclone  steering 
program . 

The  possible  modifications  of  the  numerical  forecast  procedure  according 
to  storm  stage,  path,  area,  latitude,  season,  etc.,  are  almost  limitless . 
Given  what  appears  to  be  a  suitable  numerical  steering  environment,  namely 
SR,  various  statistical  adjustments  may  now  be  derived  to  reduce  the  errors, 
especially  after  recurvature,  during  the  dissipating  stage  and  in  east-ocean 
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and  Mr.  T.  S.  Hesse,  provided  the  numerical  fields,  consulted  on  pro¬ 
gramming  and  application  and  generally  aided  in  advancing  the  research  to 
this  reporting  point.  Dr.  R.  C.  Gentry,  NHC,  Miami,  kindly  supplied 
NHC-64  data  for  comparison.  Last  but  not  least,  the  Naval  Weather  Research 
Facility,  Norfolk,  Va. ,  financially  and  morally  supported  the  project 
throughout  its  entirety. 
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TABLE  2.  -  1965  -  SR  £00  FORECAST  ERRORS  inN.autical  Miles.  Number  of  forecasts  in  parentheses 
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TABLE  3.  -  1965 

-  SR  500 
per  hour 

FORECAST  ERRORS  in 
of  forecast  interval 

nautical 

miles 

Forecast 

Interval  (hrs): 

12 

24 

36 

48 

12 

NPGS 

9.6 

8.0 

7.8 

7.3 

7.1 

Hurricanes 

OFF 

7.2 

7.1 

7.2 

7.2 

7.5 

NPGS -OFF 

2.4 

0.9 

0.6 

0.1 

-0.4 

NPGS 

6.7 

7.2 

7.1 

6.9 

6.7 

Typhoons : 

OFF 

5.8 

6.2 

6.2 

6.3 

6.1 

NPGS-OFF 

0.9 

1.0 

0.9 

0.6 

0.6 
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TABLE  4.  - 


1965  SR  700  FORECAST  ERRORS  in  nautical  miles. 
Number  of  forecasts  in  parentheses 


Hurricanes 

12  Hour 

24  Hour 

36  Hour 

48  Hour 

72  Hour 

Anna 

70(  8) 

127(  7) 

189(  6) 

284 (  5) 

524  (  3) 

Betsy 

74 (  32) 

141 (  31) 

205(  30) 

2 73 (  29) 

399(  27) 

Carol 

94(  30) 

183 (  29) 

276(  28) 

363 (  28) 

56l(  26) 

Debbie 

54  (  9) 

112(  9) 

173(  8) 

220(  7) 

261 (  5) 

Blena 

170(  14) 

328 (  13) 

471(  12) 

595 (  11) 

85  7 (  9) 

SR  700  Average 

93 (  93) 

95(  87) 

178(  89) 

183 (  83) 

262(  84) 
270(  78) 

345 (  80) 
348(  74) 

514  (  70) 

516(  63) 

SR  500  Average 

115(  89) 
114(  87) 

192(  84) 
191(  83) 

279(  79) 

2 79 (  78) 

351(  74) 

509 (  63) 

Official  Average  86(137 

171(134) 

260 

346(112) 

543 (  98) 
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TABLE  5. 


1965  SR  500  FORECAST  ERRORS  in  nautical  miles. 
Number  of  forecasts  in  parentheses 


Typhoons 

12  Hour 

24  Hour 

Lucy 

84( 

12) 

216( 

ID 

Mary 

28  ( 

7) 

UK 

6) 

Olive 

61  ( 

6) 

130( 

5) 

Rose 

34( 

9) 

139  ( 

8) 

Shirley 

122  ( 

13) 

231 C 

12) 

Trix 

97( 

14) 

187( 

13) 

Virginia 

104( 

6) 

201  ( 

5) 

Bess 

84( 

14) 

120( 

13) 

Carmen 

67( 

17) 

157< 

16) 

Della 

85( 

12) 

197( 

ID 

SR  500  Average 

80(110) 

173(100) 

80(106) 

168(  85) 

SR  700  Average 

100(117) 

180(  92) 

102(106) 

182(  85) 

Official  Average 

70 

148(199) 
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TABLE  6. - 

1965  -  12-HOUR  FORECAST  ERRORS  In  nautical  miles. 

Number  of  forecasts  in  parentheses 

Hurricanes 

Official 

SR  700  mb 

_ 

F12  F12 

Anna 

119<  6) 

o 

00 

47(  7) 

Betsy 

67(  57) 

74  (  32) 

45(  31) 

Carol 

105 (  38)  • 

94 (  30) 

42 (  28) 

Debbie 

65 (  16) 

54  (  9) 

28(  8) 

Elena 

114(  20) 

170<  14) 

71  (  13) 

Average 

86(137) 

93(  93) 

47 (  87) 

91 C  87) 


47 (  87) 


TABLE  7.  -  1965  -  12-HOUR  FORECAST  ERRORS  in  nautical  miles. 

Number  of  forecasts  in  parentheses 


SR  500  mb 

Typhoons  Official  _12 

F12  F12 


Lucy 

84( 

12) 

60( 

8) 

Mary 

28  ( 

7) 

41( 

5) 

Olive 

61( 

6) 

46( 

4) 

Rose 

34( 

9) 

43  C 

8) 

\ 

Shirley 

122  ( 

13) 

53( 

12) 

\ 

V 

Trix 

97< 

14) 

52( 

12) 

Virginia 

104( 

6) 

51( 

4) 

Bess 

84( 

14) 

54  ( 

13) 

Carmen 

67( 

17) 

62  ( 

16) 

Della 

85< 

12) 

63( 

10) 

Average 

70 

80(110) 

54(92) 

82(  92) 

54 (  92) 
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TABLE  8.  -  1965  -  24-HOUR  FORECAST  ERRORS  in  nautical  miles. 

Number  of  forecasts  in  parerttheses 


Hurricanes 

Official 

F24 

SR  700  mb 

F12 

24 

F24 

24 

Anna 

217<  4) 

127(  7) 

12  7<  6) 

126 (  5) 

Betsy 

130(  56) 

141 (  31) 

U0(  30) 

114 (  29) 

Carol 

192(  40) 

183(  29) 

108 (  27) 

1 1 7 (  26) 

Debbie 

135(  14) 

112(  9) 

72(  8) 

99<  7) 

Elena 

260(  20) 

328 (  13) 

182 (  12) 

226(  11) 

Average 

171(134) 

178(  89) 

117(  83) 

130(  78) 

174 (  83) 
173(  77) 

117(  83) 

119  (  77) 

128(  77) 
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TABLE  9.  -  1965  -  2V-H0UR  FORECAST  ERRORS  in  Nautical  Miles. 


Number 

of  Forecasts 

in  Parentheses. 

Typhoons 

Official 

f24 

SR  500mb 

pi  2 

f24 

F24 

f24 

Lucy 

169(  29) 

216(  11) 

133(  8) 

107  (  7) 

Mary 

107 (  14) 

111(  6) 

136 (  4) 

95(  3) 

Olive 

138 (  17) 

130(  5) 

5£(  3) 

124  (  1) 

Rise 

55(  15) 

139 (  8) 

147 (  7) 

58  (  6) 

Shirley 

231 (  16) 

231 (  12) 

139 (  11) 

163 (  10) 

Trix 

138 (  27) 

187  (  13) 

106(  11) 

118(  10) 

Virginia 

289  (  11) 

201 (  5) 

173  (  3) 

94  (  2) 

Bess 

106 (  26) 

120(  13) 

124 (  12) 

85(  11) 

Carmen 

148 (  21) 

157 (  16) 

134(  15) 

89 (  14) 

Della 

152 (  23) 

197 (  11) 

1 53  (  9) 

145 (  8) 

Average 

148(199) 

173(100) 

131  (  83) 

109 (  72) 

1 55 (  90) 
154 (  76) 

1 50 (  61) 

1 7 7 (  90) 

1 7 7 (  76) 
181 (  61) 

131  (  76) 

134 (  61) 

108 (  61) 

27 


TABLE 

10.-1965  -  36-HOUB 
Number 

FORECAST  ERRORS  in  nautical 
of  Forecasts  in  Parentheses 

miles. 

<  Hurricane* 

Official 

*36 

SR  700  mb 

F12 

r36 

F36 

36 

Anna 

189 (  6) 

215  <  5) 

246(  3) 

Betsy 

205 (  30) 

198 (  29) 

220(  27) 

Carpi  ^ 

276(  28) 

195 (  26) 

217(  25) 

x'  Debbie 

173<  8) 

122 (  7) 

202 (  5) 

/'  Elena 

471 C  12) 

338 (  11) 

43L(  9) 

Average 

260 

262(  84) 

211(  78) 

246(  69) 

259 (  78) 

211(  78) 

260(  68) 

215(  68) 

244 (  68) 

! 


> 

I 


f 
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TABLE  11.  -  1965  -  36-HOUR  FORECAST  ERRORS  in  Nautical  Miles 
Number  of  Forecasts  in  Parentheses 


SR  500  mb 


Typhoons  Official 

F36 

F12 

f36 

l 

F24 

36 

p12,24 

F36 

F36 

36 

Lucy 

318 (  10) 

181(  7) 

162(  6) 

152(  4) 

130(  5) 

Mary 

187(  5) 

266(  3) 

174<  3) 

226 (  2) 

152(  2) 

Olive 

mm  m 

--- 

— 

--- 

--- 

Rose 

199  (  7) 

213(  6) 

92  (  5) 

141(  5) 

100^4) 

Shirley 

335(  11) 

256(  10) 

278(  9) 

270(  9) 

308 (  8)  , 

236(  8)  v 

Trix 

275(  12) 

19  7 (  10) 

215(  9) 

180(  8) 

Virginia 

318(  4) 

341(  2) 

270(  2) 

252(  1) 

329 (  1) 

Bess  v- 

173(  12) 

194(  11) 

139 (  10) 

146(  10) 

137(  9) 

Carmen 

223(  15) 

385 (  14) 

126(  13) 

149  (  13) 

145(12) 

Della 

286(  9) 

174(  7) 

128(  6) 

104(  5) 

142 (  5) 

Average  225 

256<  85) 

247 (  70) 

170(  63) 

172(  57) 

180(54) 

255(  70) 

24  7  (  70) 

259 (  57) 

251(  57) 

162  C  57) 

259 (  57) 

251(  57) 

'  162(  57) 

172(  57) 

257(  44) 

255 (  44) 

152(  44) 

163(  44) 

169(44) 
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TABLE  12.  -  1965  -  48 -HOUR  TORECAST  ERRORS  In  nautical  miles. 

Number  of  Forecasts  in  Parentheses 

SR  700  mb  ' 

Hurricanes  Official  12  -48 

*48  *48  *48 

Anna  —  284  (  5)  304  (  4)  432(  1) 

Betsy  „ r/  273(  52)  273(  29)  313(  28)  353(  25) 

355 C  36)  363 C  28)  319(  27)  349(  24) 

310(  8)  220(  7)  198  (  6)  315(  3) 

581(  16)  595 (  11)  478(  10)  695(  7) 

Average  346(112)  345<  80)  327(  75)  391  (  60) 

341 (  75)  327(  75) 

349 (  60)  346 (  60)  39 1(  60) 
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TABLE  13. 


1965  -  48-HOUR  FORECAST  ERRORS  in  Nautical  Miles 
Number  of  Forecasts  in  Parentheses 


Typhoons 

Official 

00 

** 

F12 

F48 

SR  500  mb 

F24 

*48 

p12,24 

48 

F48 

48 

Lucy 

329 (  25) 

370(  11) 

276(  8) 

243 (  7) 

223(  5) 

267(  6) 

Mary 

154  (  10) 

293 (  4) 

364  (  3) 

190(  2) 

284  (  2) 

— 

Olive 

284(  13) 

313(  3) 

73(  I) 

— 

— 

115 (  1) 

Rose 

127(  11) 

254(  6) 

242 (  5) 

117(  4) 

173(  4) 

128 (  2) 

Shirley 

589 (  12) 

425(  10) 

383 (  9) 

379 (  8) 

364(  8) 

419 (  6) 

Trix 

'  302 (  23) 

348 (  11) 

299 (  9) 

319(  8) 

261(  7) 

438(  7) 

Virginia 

615(  7) 

418(  3) 

488(  2) 

298(  1) 

374(  1) 

— 

Beg? 

256(  22) 

229 (  11) 

278(10) 

240 (  9) 

219(  9) 

211(  7) 

Carmen 

238 (  15) 

307(  14) 

267(13) 

187(12) 

212(12) 

263(10) 

Della 

277(  17) 

366(  9) 

273(  7) 

214 (  6) 

200(  5) 

170(  4) 

Average 

304(155) 

331(  82) 

296(67) 

248(57) 

245(53) 

287(43) 

311 (  70) 

310(  59) 

314(  45) 

314(  45) 

338 (  32) 

343 (  70) 
342 (  59) 
346 (  45) 
346 (  45) 
341(  32) 

300(59) 

293(45) 

293(45) 

285(32) 

253(45) 

253(45) 

264(32) 

250(45) 

253(32) 

282(32) 
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TABLE  14.  -  1965  -  72-HOUR  FORECAST  ERRORS  in  Nautical  Miles. 

Number  of  Forecasts  in  Parentheses. 


Hurricanes 

Official 

F72 

12 

F72 

SR 

700  mb 

_24 

F72 

72 

F72 

Anna 

--- 

524 (  3) 

574  ( 

2) 

798  ( 

1) 

— 

Betsy 

469 (  48) 

399 (  27) 

564 (  26) 

563 (  25) 

553 (  21) 

Carol 

448 (  32) 

561 (  26) 

571 (  25) 

564 (  24) 

658 (  20) 

Debbie 

512 (  6) 

261 (  5) 

239  ( 

4) 

222( 

3) 

— 

Elena 

1104(  12) 

857(  9) 

888  ( 

8) 

1004  ( 

7) 

1395 (  3) 

Average 

543 (  98) 

514 (  70) 

587 (  65) 

602 (  60) 

658 (  44) 

512 (  65) 
517(  60) 

587 (  65) 

602 (  60) 

602 (  60) 
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TABLE  15.  -  1965  -  72-HOUR  FORECAST  ERRORS  in  Nautical  Miles. 

Number  of  Forecasts  in  Parentheses. 


Typhoons 

Official 

F72 

SB 

F12 

f72 

5UU  mb 

..24 

f72 

Fk 

f72 

F72 

Lucy 

481 (  19) 

526 (  9) 

332(  6) 

444 (  5) 

308 (  3) 

705  ( 

3) 

Mary 

170 (  5) 

437 (  3) 

518 (  2) 

187 (  1) 

365  (  1) 

— 

Olive 

449 (  9) 

548 (  2) 

132 (  1) 

— 

— 

— 

Rose 

245 (  3) 

453 (  4) 

431 (  3) 

148 (  2) 

268  (  2) 

— 

\ 

\ 

Shirley 

871 (  3) 

591 (  8) 

562 (  7) 

532  (  6) 

535 (  6) 

474  ( 

2) 

Trix 

426(  18) 

401 (  9) 

3  76  C  7) 

508  (  7) 

367(  6) 

706  ( 

3) 

Virginia 

1055(  3) 

786 (  2) 

690  (  1) 

— 

— 

— 

Bess 

401 (  16) 

360 (  9) 

501  (  8) 

323 (  7) 

'  360 (  7) 

265  ( 

3) 

Carmen 

289 (  11) 

441 (  12) 

466 (  11) 

329 (  10) 

397 (  10) 

516( 

6) 

Bella 

484 (  11) 

550  (  6) 

294 (  4) 

* 

237(  4) 

203  (  3) 

42  ( 

1) 

Average 

440 (  98) 

480 (  64) 

440 (  50) 

380 (  42) 

377 (  38) 

506( 

18) 

• 

458 (  50) 
463 (  38) 
463 (  38) 
483 (  14) 

440 (  50) 

446 (  38) 

446(  38) 

552 (  14) 

354 (  38) 
354 (  38) 
481 (  14) 

377 (  38) 
377 (  38) 
504 (  14) 

473 (  14) 
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TABLE  16.  -  Highlights  of  the  ova  In  it  inn  of  the 
optimal  NPGS  forecast  scheme 


i 

I 

I 


ATLANTIC: 


Forecast 

Interval 

(hrs) 

Optimal  NPGS 
Scheme 

Avg.  NPGS  Error 
(n. mi/hr  of  fest 
Interval) 

7.  NPGS 
lkt 

of  Ave. 

Errors  Within: 

2kt  3kt 

of  Ave.  of  Ave. 

Official 
Optimal  NPGS 

12 

F12 

12 

3.9 

30 

56 

76 

1.8 

24 

F12 
*  24 

4.9 

23 

46 

64 

1.4 

36 

F12 

*36 

5.9 

22 

42 

60 

1.2 

48 

F12 

*48 

6.8 

18 

40 

55 

1.1 

72 

F72 

7.1 

17 

34 

50 

1.1 

PACIFIC: 

12 

F12 

*12 

4.5 

28 

60 

82 

1.3 

24 

F24 

24 

4.5 

34 

51 

74 

1.5 

36 

F24 

*36 

4.7 

30 

54 

69 

1.2 

48 

F12,24 

48 

5.1 

21 

46 

64 

1.2 

72 

v12'24 

*72 

5.2 

21 

41 

62 

1.2 

Error 
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TABLE  18.  -  1965  ATLANTIC  HURRICANES:  ERROR  DISTRIBUTION  FOR  OPTIMAL  NPGS  SR  700  SCHEME  BY  AREA,  PATH  AND  STORM  STAGE 

Errors  In  Nautical  Miles.  Number  of  Forecasts  in  Parentheses 


/*> 

X*N 

x-s 

x-s 

CM 

CM 

CO 

CO 

•W 

K-/ 

S/ 

w 

s-x 

H 

O 

m 

r*l 

ON 

O 

O 

CM 

m 

m 

r* 

IB 

r-l 

CM 

CO 

SO 

SI 


no 

cn 


x-s 

xn 

X** 

x-s 

vO 

vO 

vO 

sO 

s_x 

S.X 

v-x 

s-x 

00 

CO 

Os 

Os 

CM 

Os 

so 

H 

m 

00 

XN 

X*\ 

H 

CM 

rM 

o 

m 

m 

m 

VO 

M/ 

s-x 

V 

M/ 

vO 

so 

o 

r-l 

r-l 

CM 

00 

rM 

CM 

cn 

St 

XN 

XS 

x-s 

x-s 

CO 

H 

ON 

r^ 

f— 1 

s-x 

sx 

S-x 

SX 

§ 

Os 

O 

00 

cn 

3 

r-l 

CM 

CM 

<t 

xs 

X-N 

X-N 

x-s 

in 

r-l 

oo 

vO 

<t 

r-< 

H 

sx 

V-X 

nx 

V-/ 

sx 

n* 

CM 

CM 

On 

cn 

H 

H 

H 

CM 

st 

x-v 

x-s 

<t 

in 

St 

m 

St 

CM 

CM 

CM 

CM 

CM 

S-x 

s-x 

w 

v-< 

\x 

ON 

<t 

O 

O 

so 

so 

<JV 

co 

00 

r-l 

CM 

vO 

a) 

X-s 

x-\ 

x-s 

XN 

XN 

n 

CO 

00 

o 

VO 

o 

sO 

in 

in 

in 

st 

MS 

s-x 

S-x 

s-x 

NX 

S-X 

ON 

m 

vO 

IN¬ 

£ 

CO 

ON 

rs 

IN 

H 

CM 

st 

X"N 

x-s 

X-N 

XN 

XN 

1 f 

m 

m 

m 

in 

U  '  1 
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s-x 

sx 

NX 

O 

SO 

Ov 

m 

's. 

CO 

in 

CM 

rs 

H 

CM 

M0 

X-N 

x-s 

X-N 

XN 

x-s 

St 

CO 

r-^ 

ON 

r-» 

CO 

CO 

CO 

CM 

CM 

nx 

S-x 

S-x 

NX 

sx 

CM 

r^ 

Ctv 

cn 

St 

os 

r^ 

vO 

CM 

rH 

CM 

cn 

x-s 

x-s 

x-s 

XN 

XN 

00 

in 

CM 

•■M 

00 

St 

cn 

<t\ 

s-x 

s-x 

S-X 

MX 

s^ 

O 

H 

m 

NO 

00 

m 

CO 

CO 

m 

CM 
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CM 

<n 

M3 
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CM 
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TABLE  15.  -  1965  PACIFIC  TYPHOONS:  ERROR  DISTRIBUTION  FOR  OPTIMAL  NPGS  SR  500  SCHEME  BY  PATH 

Errors  In  Nautical  Miles.  Number  of  Forecasts  In  Parentheses 


CM  CM 

?-<  p* 

Em  • 


pm" 


<t  vO 

cm  n 
Pm 


•Jt  00 

%m* 


CM 
CM  f" 
PM 
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CAPTIONS  FOR  FIGURES 


FIGURE  la  -  Portion  of  FNWF's  operational  500-mb  height  (Z)  and  temperature 
analysis  for  0000  GMT  20  August  1965.  Contours  at  60-meter  interval 
(dark  solid  lines);  isoline  (center)  labels  in  tens  (units)  of  meters  with 
thousands  figure  omitted.  Isotherms  (light solid  lines)  not  labeled. 

FIGURE  lb  -  Portion  of  FNWF's  operational  500-mb  analysis  of  the  disturb¬ 
ance  component,  ZgD,  derived  from  the  500-mb  height  field,  Z,  at 
0000  GMT  20  August  1965.  Contours  at  30-meter  interval;  isoline 
(center)  labels  in  tens  (units)  of  meters . 

FIGURE  lc  -  Portion  of  FNWF's  operational  500-mb  analysis  of  the  residual 
or  smoothed  component,  ZgR,  derived  from  the  500-mb  height  field, 

Z,  at  0000  GMT  20  August  1965.  Contours  at  60-meter  interval;  labels 


as  in  figure  la . 

FIGURE  2  -  Schematic  representation  of  Z 


SR 


contours  with  superimposed 


FNWF  linear  grid.  I,  J  is  tropical  cyclone  location.  Ge9strophic  winds 
(V  )  ,are  computed,  by  component  (V  ,  V  ),  at  the  four  identified  grid 
points,  each  at  a  distance  D  from  I,  J. 

FIGURE  3  -  Graph  of  sin  0  and  mod  sin  0  =  2[(.25  sin  «  +  .25)  +  .25  sin  0  ]. 
FIGURE  4  -  Hurricane  Elena  positions  at  12-hour  intervals,  starting  at 

1800  GMT  12  October  1965  (position  "0"):  best  track  (•* - *)  and 

numerical-steering  forecast  track,  F24  (—7-—  ). 

FIGURE  5  -  Schematic  example  of  a  modified  24-hour  numerical  steering 

forecast  (using  E^^)  made  from  best  track  position  T  .  In  general,  the 

—  ^ 

vector  correction  for  the  bias  in  numerical  steering,  E  ,  is  for  the 

yy 

same  time  period  as  the  forecast  interval  implied  by  F  and  results 
in  the  modified  position  designated  as  F  ,  where  yy  =  xx. 

FIGURE  6  -  Schematic  example  of  a  modified  24-hour  numerical  steering 
forecast  (using  E^)  made  from  best  track  position  T  .  In  general, 
the  vector  correction  for  bias  in  numerical  steering  is 
forecast  interval 


E12x 


12 


and  results  in  the  modified  position  designated 


as  F 


12 

xx 
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FIGURE  7  -  Same  as  figure  4  with  the  addition  of  forecast  tracks 
(  —  •  —  A  —  •  — )  and  (  •  •  0  •  •  )  . 

FIGURE  8  -  Schematic  example  of  a  modified  36-hour  numerical  steering 

In  general, 

forecast  interval 


forecast  (using  E„  .}  made  from  best  track  position  T 
Z4  o 


the  vector  correction  for  bias  in  numerical  steering  is  E 24  x 


24 


and  results  in  the  modified  position  designated  as  F 


24 

xx 


FIGURE  9  -  Schematic  example  of  a  modified  36-hour  numerical  steering  fore- 

cast  (using  E,  „  and  E„  .)  made  from  best  track  position  T  .  In  general, 
1Z  Z4  o  . 

the  vector  correction  for  bias  in  numerical  steering  is  a"?12  +  BE^  where 


,  forecast  interval 

A  =  - — -  and  B 

24 


.5A,  and  results  in  the  modified  position 


_12,24 

designated  as  F 

XX 

FIGURE  10  -  Typhoon  Carmen  positions  at  12-hour  intervals,  starting  with 

0600  GMT  1  October  1965  (position  "0"):  best  track  (-* — if-)  and 

48 

modified  numerical-steering  forecast  tracks,  F4g  (•  —  A  —  *), 

(  —  V  — ),  andF*2'24  (..0*.)  . 

FIGURE  11  -  Tropical  storm  Debbie.  12,  24,  36,  48  and  72  hour  forecasts 
made  at  1800  GMT  25  September  1965,  using  the  numerical- steering 
computation  F  (  —  7  — ),  and  optimal  scheme  Fopt  (•  •  0  *  • )  . 

XX  XX 

Best  track  (-#* — »t-)and  available  official-forecast  positions  (  A) 
are  shown. 

FIGURE  12  -  Typhoon  Rose,  forecasts  made  at  1800  GMT  2  September  1965. 
Remainder  of  legend  as  in  figure  11. 

FIGURE  13  -  Division  of  North  Atlantic  area  (A,  B,  C)  used  in  stratifying 
forecast  statistics  in  table  18. 
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